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Abstract 

For a paradigmatic case, the standard map, we discuss how the statisti- 
cal description of the approach to equilibrium is related to the sensitivity 
to the initial conditions of the system. Using a numerical analysis we 
present an anomalous scenario that may give some insight on the founda- 
tions of the Tsallis' statistical mechanics. 
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1 Introduction 

Since long time it is a widely diffused opinion between physicists that the basis 
of statistical mechanics lies on dynamics. Particularly, Krylov has pointed 
out that the mixing properties of a dynamical system are responsible for its 
statistical behaviour (herein we use this word as a synonym of sensitivity to 
initial conditions: i{t) = limAa:(o)->o ^2;(t)/Aa;(0), for a one-dimensional illus- 
tration). Along this line, intensive work has been done, especially in situations 
were anomalous effects may arise (see, for example, Q and references therein). 
In it was studied numerically the approach to equilibrium of an Hamiltonian 
system, the standard map, also referred to as the kicked-rotator model: 



xt+i = Vt + sin(27rxt) -I- Xt (mod 1), 

ZTT 



Vt+i = + 7^ sin(27ra;t) (modi), 

ZTT 



(1) 



where a 6 M. This map presents a simplectic structure that corresponds to an 
integrable system when a = 0, while, for large-enough values of a, it is strongly 
chaotic (in Fig. ^ we display the phase portrait of the map (^ for typical values 
of a) . It was shown that for small values of the parameter a (where the border 
between the regular and the chaotic region becomes significant), the approach 
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Figure 1: Phase portrait of tlie standard map for typical values of a. Af = 20 x 20 
orbits (black dots) were started with a uniform ditribution and traced for < t < 200. 



to equilibrium of the system displays an anomalous behaviour of the usual, 
Boltzmann-Gibbs (BG), statistical entropy; anomaly that may open the door 
for a dynamical foundation of the Tsallis' thermodynamics. In this paper we 
review the numerical analysis |^ , presenting some new results and making some 
further speculations. 



2 Statistical description of the approach to equi- 
Ubrium 

Assuming a Gibbsian point of view, we can study the system (|^) (in its F space) 
proceeding, for example, as follows. We start introducing a coarse-graining 
partition of the phase space by dividing it in W cells of equal size, and we 
set many copies of the system {N points) in a far-from-equilibrium situation 
putting all the TV points inside a single cell. The occupation number Ni of each 
cell i — ^) provides a probability distribution pi — Ni/N, hence the 

definition of an entropy value: 

-. q 

We remind that the entropic form reduces to the BG entropy Si = — Pi In Pi 
in the limit q ^ 1 (for a recent review on Tsallis' statistics, see [Q). Using then 
the dynamic equations (^, at each step the points spread in the phase space 
causing the entropy value to change. Finally, to extract a global quantity on the 
mapping phase space, we repeat the calculation setting the cell that contains 
the initial points in different positions chosen randomly all over the whole unit 
square, and we take an average over all the different histories thus obtained. 
The key point is to perform many different histories so that the average stabi- 
lizes on a definite curve. The result of this analysis for fixed a is then a single 
curve of the entropy versus time, for each entropic form 5*^ (in Fig. ^ we show 
a typical case). 



2 




Figure 2: Sq{t) for a = 2 and q = 1, with N = W = 1000 x 1000 and initial data 
chosen randomly all over the whole unit square (see text), (a) 50 different histories, 
(b) average over 5000 histories. 

3 Power-law mixing as a foundation of Tsallis' 
generalized statistical mechanics 

We know from the theory of chaos that when a map displays a strong chaotic 
behavior, the mixing is exponential. When we attempt an analysis like the 
one described in the previous section, there is just one value of q for which the 
entropy displays a linear stage before saturation: q = 1; moreover (see [Q), the 
slope of this linear stage is equal to the Kolmogorov-Sinai entropy rate (i.e., 
the positive Lyapunov exponent). The resultant curve has the features of Fig. 
^ (b). If in the same time interval where Si(t) grows linearly we use q < I 
{q > 1) the curve bends upward (downward) before saturation. This situation 
recalls the definition of the Hausdorff's dimension of a space: taking first the 
Umit W ^ oo and then the limit i — > cx) there is only one value of q for which 
the entropy production is different from and from +oo. 

On the other side, if the phase space has characteristics like those in Fig. 
|l|(c), the islands-around-islands structure at the border between the strongly 
chaotic and the regular regions (see, for example, Q) slows down the mixing, 
making it to become a power-law mixing instead of an exponential one. In 
this case, for time not too large, the linear growth with time occurs for the 
entropy Sg* with g* < 1 {q* ~ 0.3 for the standard map), and not for 
(see Fig. y(a) and ^(c)). Waiting enough time, after t = tcross, a crossover 
to the exponential mixing would occur, due to the rapidity of the exponential 
growth with respect to the power-law one. We have then two different statistical 
regimes for < \a\ < 1: first (t << tcross)i one for which the mixing properties 
are well described by a Tsallis' entropy Sq<- with q* < 1; then {t >> tcross), one 
for which the mixing properities are well described by the usual BG entropy. 
The interesting point is that for some particular macroscopic conditions (here 
represented by the value of the parameter a) , the complexification of the phase 
space (especially when dimensionality increases) may increase tcross, keeping 
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Figure 3: Sq{t) for typical values of a (5000 to 16000 histories were averaged); squares 
correspond to N = W = 2236 x 2236; circles correspond to iV = = 1000 x 1000. 
The lines are guides to the eye. (a) g = 1, (b) q = 0.5, (c) q = 0.3, (d) q — 0.05 
The arrows indicate the crossover time tcross for the corresponding values of a, as 
represented in the inset of (c). 



the system, for at least some interesting physical observations, in a meta-stable 
Tsallis' phase. This scenario is consistent with some anomalies observed in long 
range many-body Hamiltonians |^ that were presented at this congress. 

In Fig. H we present, for typical values of q, how this anomalous stage in 
the entropy growth makes its appeareance, when one decreases \a\ under certain 
value. The crossover time tcross increases when a — > (see inset of Fig. ||(c)). 
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